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ABSTRACT 

We present spectroscopic rotation velocities (v sin i) for 56 M dwarf stars using high reso- 
lution HET HRS red spectroscopy. In addition we have also determined photometric effective 
temperatures, masses and metallicities ([Fe/H]) for some stars observed here and in the litera- 
ture where we could acquire accurate parallax measurements and relevant photometry. We have 
increased the number of known v sin is for mid M stars by around 80% and can confirm a weakly 
increasing rotation velocity with decreasing effective temperature. Our sample of v sin is peak 
at low velocities (~3 km s^ 1 ). We find a change in the rotational velocity distribution between 
early M and late M stars, which is likely due to the changing field topology between partially 
and fully convective stars. There is also a possible further change in the rotational distribution 
towards the late M dwarfs where dust begins to play a role in the stellar atmospheres. We also 
link v sin i to age and show how it can be used to provide mid-M star age limits. 

When all literature velocities for M dwarfs are added to our sample there are 198 with 
v sin % < 10 km s _1 and 124 in the mid-to-late M star regime (M3.0-M9.5) where measuring 
precision optical radial- velocities is difficult. In addition we also search the spectra for any sig- 
nificant Ha emission or absorption. 43% were found to exhibit such emission and could represent 
young, active objects with high levels of radial- velocity noise. We acquired two epochs of spectra 
for the star GJ1253 spread by almost one month and the Ha profile changed from showing no 
clear signs of emission, to exhibiting a clear emission peak. Four stars in our sample appear 
to be low-mass binaries (GJ1080, GJ3129, G1802 and LHS3080), with both GJ3129 and G1802 
exhibiting double Ha emission features. The tables presented here will aid any future M star 
planet search target selection to extract stars that will exhibit low radial-velocity jitter. 

Subject headings: stars: fundamental parameters — stars: low-mass, brown dwarfs — stars: rotation — 
(stars:) planetary systems 



1. Introduction 

In the past 10 years radial-velocity measure- 
ments in the optical have made great strides by 
utilising a number of techniques and methodolo- 



gies to gen erate precisions of 3 m s _1 in the long 

and sub-m s _1 in 
Use of 



term (e.g. iButler et al.l 120061 



the short term (e.g^ 

an iodine cell ( e.g. iMarcv fc Butlerlll992Tl or ThAr 



Bouchv et al 



2005). 



gas lamp (e.g. IPepe et ah 2000l ) have allowed de- 



1 Based on observations obtained with the Hobby-Eberly 
Telescope, which is a joint project of the University of 
Texas at Austin, the Pennsylvania State University, Stan- 
ford University, Ludwig-Maximilians-Universitat Miinchen, 
and Georg-August-Universitat Gottingen. 



tection of around 300 extrasolar planets (exoplan- 
ets), with that number increasing each month (see 
http://exoplanet.eu/). A large parameter space 
has been studied but since these observations are 
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limited to the optical regime, where M dwarfs 
are intrinsically faint, a vast amount of stars are 
left unobserved, particularly the M star popula- 
tion that constitutes the bulk of stars in the local 
galactic neighbourhood. Radial-velocity studies of 
early M dwarfs have detected planets well into the 
terrestrial-mass regime, down below 1 0M m (e.g . 



reg 

Rivera et alJ 12003 : lUdrv et alJl2007t iMavor et al 
2009h . There is considerably improved mass con- 



trast with such stars (amplitude of a given mass 
planet oc M s t a r ' 5 ) and obtaining precision radial- 
velocities of these objects may enable Earth-mass 
planets to be detected in their habitable zones. 

The largest uncertainty associated with preci- 
sion radial-velocity measurements is that of stel- 
lar activity. The association between stellar ac- 
tivity and rotation velocity (v sin i) is well es- 
tablished ( Noves et al. 1984h and measuring this 
parameter is an excellent proxy of the level of ac- 
tivity. For a fixed resolution, S/N and calibra- 
tion method, stars with higher rotational velocities 
have radial- velocity measurements with lower pre- 
cision since the larger rotation serves to w ash out 
the spectral features (jBouchv et al.1 120011 ). Cou- 
pled to this, the combination of activity and rota- 
tion can cause false pos itives to appear in the data 
(e.g. I Henry et al. I l2002h but in general it s erves to 
increase the level of jitter (Wright] 2003Y How- 
ever, the rotation-activity connection has been 
shown to saturate quickly aroun d M stars, occur- 
ingat ~5 km s _1 for M2 objects ( Patten fc Simonl 
1996). Since we are probing later M stars (>M3), 
where observations of various activity indicators 
have shown that the percentage of active stars 
will increase drama t ically, up to ~100% at M7 
(jFleming et alj|200d : iMohantv et~ai1l2002f) , we are 
focused on selecting the slowest rotators within 
each spectral bin, in order to select the narrowest 
line profiles. An increase in radial- velocity uncer- 
tainty of around 3-4 times was found by Bouchy 
et al. when increasing the v sin i of stellar models 
from 1-10 km s _1 . The larger uncertainty arises 
due to the loss of spectral information in the stars 
with higher v sin i since the blending factor is in- 
creased within the stellar forrest, to the detremcnt 
of the radial- velocity information. 

Based on the above description any future pre- 
cision radial-velocity planet search survey target- 
ing cool M stars should have a fairly strict selec- 
tion based on the rotation velocity of their sam- 



ple stars, particularly when targeting Earth-mass 
planets in the habitable zones. For instance, a 
search employed using a PRVS-like instrument 
(| Jones et al J 120081; iRamsev et al1l2008l) would fo- 
cus on selecting the brightest M stars with ro- 
tational velocities in the sub-10 km s _1 regime, 
which would allow high S/N spectra to be ac- 
quired in the shortest possible observational times, 
gaining in the number of M stars observed in a 
single observing run, and also allowing the high- 
est precisions to be reached to detect the lowest 
mass rocky planets, in particular Earth-like plan- 
ets in their habitable zones. We do note that 
such a 10 km s^ 1 upper limit should probably 
be raised when moving into the ultracool/brown 
dwarf regime to increase the number of such ob- 
jects on any planet search sample since previous 
studies have shown that rotatio n velocities appear 
to systematically increase ( e.g. Mohanty k. Basril 
20031 ; iReiners fc BasrilfeoOSl ). 



2. Observations & Reduction 

All observations in this study were made over 
the period 2006 December 25-2007 November 04 
utilising the queue scheduling mo de at the 9.2m 
Hobb y-Eberly Telescope (HET; IRamsev et al 



1998) at McDonald Observatory in Texas. Us- 
ing t he High Resolution Spectrog raph (HRS;[TuI3 
19981 ) red chip (~8100-9900A) operating at a res- 
olution of i?~37,000 and spatial binning of 2x2, 
53 objects were observed, ranging in spectral type 
from M3V to M6.5V. The HRS observing mode 
employed the use of two 3" fibers since the objects 
were so faint, with one on object and one on sky. 
Bright telluric standards taken from the HET list 
of rapidly rotating B stars were also observed on 
most of the nights, which allowed us to charac- 
terize the level of telluric contamination in the 
data. Since even close by mid-to-late M stars are 
relatively faint, our HRS integration times ranged 
in duration from 5-25 minutes and for some of 
the fainter objects multiple exposures were taken 
and combined to generate a high S/N spectrum. 
However, most of the v sin i measurements, when 
comparing the individual frames to the combined 
frame, were within ±0.5 km s _1 since the decon- 
volution method we employ boosts the S /N in the 
final line profile used in the measurements. 

The reduction of all data was performed using 
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CCDPACK and STARLINK techniques. First the bias 
and overscan signals were removed and the over- 
scan region was clipped off each image. A master 
ffatfield was created by median filtering all flats 
taken in the standard HRS calibration plan. This 
usually consisted of a few bias frames, flatficlds 
and ThAr arc observations taken before and after 
the nights observing, however for one night no flats 
were observed and a flatfield from a night close 
to the observing night was used in the reduction. 
Also the number of flatfields were increased to 10 
for the last six months or so of data, which signifi- 
cantly aids in fringe removal. Indeed, the fringing 
in the red chip of HRS can reach 10-20% and hence 
well exposed, high S /N Halogen flatfields are used 
to correct these. 

Th e starlink package echomop ([Mills et alj 
1996) was used to extract the echelle data. Firstly, 
all 14 orders were located and traced using the 
master flatfield, and any stray pixels were clipped 
from this trace. The dekker limits were deter- 
mined using the master flatfield to ensure the en- 
tire object order and the adjoining sky order are 
included in the reduction. The master flatfield was 
then normalized to determine the balance factors 
and these were multiplied into all other observa- 
tions. The sky light was removed by utilising the 
sky fiber and subtracting out any sky lines that 
appeared in the object spectra. The orders were 
then extracte d using an optimal extraction routine 
( Home Il986l) and wavelength calibrated using the 
ThAr arc images. 

A small portion of spectra for the M4.5V star 
G 121-028, centered on the sodium doublet at 
8182-8194A, is shown in Fig. CD (solid line). The 
dashed spectrum is the telluric standard HD1839, 
which is a fast rotating B star, therefore all the 
strong features are due to telluric contamination. 
The sodium doublet is heavily blended with a 
number of telluric features, which without good 
sky subtraction, could serve to contaminate any 
measurement of the line properties. This is one 
of the reasons that we use a spectral deconvolu- 
tion routine to help remove the telluric features 
in our final analysis. These final spectra have a 
limiting resolution of between ~0.20-0.35A (~7- 
13 km s _1 ), which was determined individually 
for each star by deconvolving the telluric lines in 
each image. 



3. v sin i Determination 

All v sin is in this work are determined by deriv- 
ing an optimal line profile for each star, and then 
convolving a non-rotating template (LHS1950 is 
our template) with rotational profiles of various 
velocities to find the best match. A deconvo- 



lution method (IDonati & Collier Cameron 


19971: 


Barnes et al. 1998 


; Collier Cameron et al. 


2002) 



is used to determine the optimal line profile. A 
continuum fit is performed using the CONTINUUM 
routine in iraf to normalize the spectra and to 
provide inverse variance weights for the deconvolu- 
tion procedure. This least squares deconvolution 
uses a line list which best matches the spectral 
type of each star. The wavelength positions and 
line depths used were take n from VALD (Vienna 



Atom ic Line Database; iKupka fc Rvabchikova 
1999). Any deficiencies in the line list data for 
cool M stars, such as missing opacities, should not 
affect the final v sin i values since the same list is 
used on the non-rotating template star LHS1950. 
The use of fewer model lines for deconvolution 
than actually observed will simply give rise to 
profiles with lower S/N than a more complete 
model would afford. The line strengths of the fi- 
nal profiles may also vary from star to star but 
can be corrected for by scaling the non-rotating 
template to the line strength of each stellar de- 
convolved profile. The method of least squares 
deconvolution is used to determine the mean line 
profile which when convolved with the line depth 
pattern gives the optimal match to the observed 
spectrum. 

Fig. [2] represents the final deconvolved line pro- 
file for G 121-028 (solid line) with velocity on the 
x-axis and normalized flux on the y-axis. G 121- 
028 is one of the slower rotating stars in our sam- 
ple (v sin i=3.8±0.7 km s _1 ) as the narrow line 
profile shows. To determine the v sin i values 
we perform a Levenberg-Marquardt l east squares 
minim isation using MPFIT in IDL (jMarkwardt 
2009) to find the best fit for each observed line 
profile. The fit is represented by the dashed line 
in the figure and for this example a reduced \ 2 
(Xv 2 ) of 1-2 is found. This best fit is used to de- 
termine both the profile centroid and scaling fac- 
tor. The centroid allows each star to be shifted to 
the rest frame for use in the template comparison, 
whilst the scaling factor allows us to scale our non- 
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Fig. 1. — The normalized HRS spectra for the 
M4.5 dwarf G 121-028 around the sodium dou- 
blet (8182-8194A) is shown by the solid line. 
The dashed line represents the telluric standard 
HD1839 and has been offset for clarity. Note the 
blending between the strong sodium lines and the 
tellurics in the M dwarf spectra. 
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Fig. 2. — The deconvolved line profile for the star 
G 121-028 in velocity space (vacuum velocities). 
The best fit to this profile is marked by the dashed 
line. 



rotating template to match the profile strength of 
each star. LHS1950 was found to be our narrow- 
est profile (determined from the full width at half 
maximum (FWHM) of the best fit to each profile) 
and hence below the resolution limit of the instru- 
ment. This template is then broadened in steps of 
0.5 km s _1 , between velocities of 0.5-50 km s" 1 , 
using a grid of rotational profiles with a limb dark- 
ening coefficient of 0.6. An example is shown in 
Fig. [3] A zoomed in region of the profile of G 121- 
028 is represented by the dotted curve, with tem- 
plate profiles marked by the solid curves. From 
inner to outer (narrowest to broadest) the solid 
curves have been broadened by rotational profiles 
of 3.5, 5.0 and 10.0 km s _1 respectively. This 
highlights the good fit between the inner profile 
compared to the outer profiles within the uncer- 
tainties and highlights the robust nature of the 
fitting technique. All rotational profiles are g ener- 
ated following the description in lGravl (|1992l ). and 
a \ 2 value is determined at each v sin i step. 

Fig. [4] shows our x 2 fit to the data points in 
the analysis of the star G 180-011. Each filled cir- 
cle represents one of the 0.5 km s^ 1 steps used to 
broaden the template to fit the profile of G 180- 
011. The solid curve represents the best quadratic 
fit to the data, an d following the procedure in 
Jenkins et al. ( 2008h . the minimum of this curve is 



the measured v sin i, highlighted here by the verti- 
cal solid line. The vertical dotted, dot-dashed and 
dashed lines represent the la, 2a and 3tr inter- 
nal uncertainties respectively. The la uncertainty 
value was determined by measuring the width of 



the curve at a la step (i.e. 



X 2 ia 



^ 7/ 



+ 1), 



after the curve has been broadened by a factor 
determined from the difference between the mini- 
mum of the x„ 2 and 1 (where \ v 2 is the best fit). 
Due to the nature of the dataset this factor has a 
median of only 2.1 and standard deviation of 1.5, 
indicating the goodness of the methodology. We 
employ this artificial broadening to the x 2 curve 
to help alleviate uncertainties that are difficult to 
address in the analysis procedure e.g. macrotur- 
bulence variations, spectral type differences, etc. 
For G 180-011 shown in the figure the formal un- 
certainty is found to be ±0.4 km s _1 . We use this 
procedure to determine all la uncertainties quoted 
in Table [3] 

The la uncertainties mentioned show we have 
good internal precision and since we have a few 
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Fig. 3. — A zoomed in region of deconvolved pro- 
file of G 121-028 (dotted curve), along with three 
broadened template profiles (solid curves). From 
inner to outer the template has been broadened by 
velocities of 3.5, 5.0 and 10.0 km s -1 respectively. 



stars where multiple observations were taken we 
can test this in a brute force manner by looking at 
the values for these multiple measurements indi- 
vidually. When we do this we find that the v sin is 
tend to agree to within the uncertainties deter- 
mined by the y 2 fitting procedure on each indi- 
vidual measurement. This gives us confidence in 
the analysis procedure. We note that for the star 
GJ1253 the v sin i values agree to within ler, which 
for this star was ~1 km s _1 , and as these were 
measured over two observations a month apart, 
they are probably the best indicator of the over- 
all random uncertainties in the analysis procedure. 
Therefore, all velocities for objects <10 km s _1 
are accurate to ±1 km s _1 . For the objects ro- 
tating much faster than this, particularly above 
20 km s _1 , it is difficult to fit their profiles to bet- 
ter than ~±10-20% at the adopted resolution and 
S/N and therefore this should be taken as the ac- 
curacy for such fast rotators. 
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Fig. 4. — A subsection of the best fit \ 2 curve to 
the data (filled circles) for the star G 180-011. The 
solid curve marks the best fit, with the minimum 
of this function marked by the vertical solid line. 
The observed rotational velocity is 21.9 km s _1 
for this star. The dotted, dot-dashed and dashed 
lines mark the determined la, 2er and 3<r limits 
respectively. 



3.1. Instrumental Profile 

Although we estimate our resolution limit to be 
around 2.5 km s _1 , we note that this is not the case 
for all of our sample. Since we use a fixed template 
spectrum that was measured on our first night and 
not on every night of our observations, there will 
be effects due to the changes of the instrumental 
response. To monitor this we used the telluric lines 
in each of our spectra as a proxy for the instru- 
mental profile (IP). We assume that the intrinsic 
width of weak telluric lines are smaller than the 
instrumental resolution. We then performed the 
same fitting routine to these deconvolved profiles 
which returned an instrumental FWHM for each 
stellar spectrum. We decided not to use the arc 
lines for the IP measurements since the calibra- 
tions are taken through a different optical path 
than the science images and therefore would not 
acc urately mo del the intrinsic instrumental width 
(see lTulll998l for HRS design). This procedure al- 
lowed us to determine that there was an increase 
in the instrumental width of almost a factor two 
between 2007 August 9 th and 2007 October 17 th 
due to poor thermal control in the HET spectro- 
graph that lead to PSF instability through focus 
drift. The stars observed between these dates have 
a typical resolution limit of around 4-5 km s _1 , 
therefore we employ a correction to each star to 
correct for the effects of the broadening of the IP. 
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Fig. [5] shows the average change in the FWHM 
of a selection six stars with the highest S/N ra- 
tios as a function of v sin i. Each of the data 
points mark the average of the FWHM difference 
(FWHM-FWHMo) for the six stars, along with the 
associated scatter of the values, after broadening 
by rotation profiles in steps of 0.5 km s _1 , where 
FWHM is the initial FWHM before any broaden- 
ing has been applied. The solid line shows the best 
straight line fit to the data that we employed to 
correct for the IP, and is described by 1.08x-2.03, 
with an RMS scatter of only 0.46 km s _1 . 

To employ these corrections we determine the 
difference between the template telluric FWHM 
and that of each star. We then use the fit to 
determine which v sin i profile we must broaden 
the template profile by initially to correct for the 
changing IP width. Due to our resolution this 
only becomes significant on the data mentioned 
above where the IP almost doubled in width, leav- 
ing such stars with a resolution of 4.5 km s _1 . 
These fits can also be used to determine the actual 
v sin i from the change in the FWHM compared 
to the template, at least in the low v sin i regime, 
and a test of this reveals agreement between both 
methods. The fits could be extended into the high 
v sin i regime in order to use this method to de- 
termine the rotation velocities for rapidly rotating 
stars also. All v sin is and their associated uncer- 
tainties are shown in column 1 1 of Table [3l with 
the telluric FWHMs shown in column 12 for ref- 
erence. 

Fig. [6] shows the overall spread in the measured 
v sin i values as a function of the changing IP of 
the instrument. As mentioned above the change 
in the IP is monitored by measuring the FWHM 
of the telluric lines in each spectrum (FWHMxeii) 
and from the plot there appears no significant cor- 
relation between measured rotation values and the 
width of the IP. The linear correlation coefficient 
(r) for the entire sample is -0.12 and is still only 
0.36 when removing all stars with measured upper 
limits, which are represented in the figure by the 
downward pointing arrows. This test highlights 
the lack of any significant correlation between the 
measured v sin is and the changing width of the 
IP. In fact, only one of the stars measured dur- 
ing the period of increased IP values is found to 
be above the detection threshold of our method, 
and as will be seen later, this correlates with the 




Fig. 5. — The change in the FWHM as a function 
of v sin i. The squares show the average FWHM 
minus the initial FWHM (FWHM ) after the pro- 
file has been broadened in steps of 0.5 km s . 
The solid curve is the best straight line fit to the 
data. 
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Fig. 6. — The distribution of v sin i as a function 
of the change in width of the instrumental profile. 
The width is parameterised by the FWHM of the 
telluric lines in km s _1 . No correlations are evi- 
dent in the spread meaning the correction applied 
to the template spectrum appears robust. The 
symbols with downward pointing arrows indicate 
upper limits. 
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presence of Ha in emission in this star. 

None of these objects have previously deter- 
mined v sin is with which to c ompare our values 
for an accuracy check. However. iRockenfeller et al. 
(|2006l ) determined the /-band variability for the 
star LHS2930 and found the period to be 13±2hrs. 
They also quote the radius of the star to be 
O.33R,0, from which we can estimate the rotation 
velocity assuming the variability is induced by the 
rotation of surface features such as star spots. Es- 
timating this velocity in pure spherical geometry 
returns a value for the equatorial rotation veloc- 
ity of S0t 6 3 km s _1 . We measure a velocity of 
18.7±1.5 km s"\ which is significantly lower than 
the photometrically derived value. We note that 
the deconvolved line profile for this star is flagged 
as low S/N. However, given these values this would 
indicate the star is inclin ed to our line of s i ght b 



around 51I.!- Uo («~39°). iMarcv fe Butler! (|199 



show that for random stellar alignments, 68% of 
stars will have inclinations to our line of sight of 
below 47°, with 95% having inclinations above 
~20°, placing LHS2930 in between this 1 and 2a 
result and in agreement with the la value within 
the estimated uncertainty. 

3.2. Ha Absorption/Emission 

In addition to selecting against stars with a 
large v sin i, one might also like to monitor each 
star's magnetic activity through the absorption or 
emission in the Ha line. The correlation between 
Ha emission, age, activit y and rotation velocity 
h as been well studied ( e.g. 



/ ana rot ation velocity 
Mohantv fc Basril2003t 



Reiners fc Basril 120081 and refs therein), however 
the magnetic field structure of M dwarfs, partic- 
ularly slowly rotating M dwarfs, could be com- 
plex, with magnetically active regions rotating in 
and out of view and the presence of strong flar- 
ing events. Also there is the possibility that a 
small percentage of our M dwarfs are being viewed 
close to pole-on and instead of measuring their 
true v sin i, which could be very high, we measure 
a smaller value and place these in the good planet 
search target bracket. Therefore, we searched for 
any emission, and also significant absorption, in 
the Ha line (A~6562A) in each of our target stars. 
24 (43%) stars were found to show Ha emission, 
which includes the double profiles found for the 
active binary stars, and since such stars could 
be young and highly active, these may not be 



ideal stars to include in a precision radial- velocity 
planet search. Of these only six have rotation ve- 
locities below 10 km s _1 (12% of the total). How- 
ever, we do caution that we might be viewing ac- 
tive regions on these stars, or flaring outbursts, 
since numerous M stars across the spectral do- 
main have b een found to exhibit such phenom- 
ena (e.g. see Tinncv ct al.lll998t Martin fc Ardila 
200l|; lOsten et all 20051 ). Further measurements 



may reveal these to vanish and therefore more 
epochs may be required to test if such stars have 
continuous Ha emission and are probably highly 
inclined young and active sta rs, or are simpl y go- 
ing through a flaring event. Reinerd ( 20091 ) has 
shown that flares only create significant noise at 
velocity precisions below 10 m s _1 for moderate 
events and at the level of a few hundreds of m s _1 
for giant flares. However, giant flares are also 
heavily correlated with Ha variation and so can 
be easily removed from a radial- velocity campaign, 
albeit at the cost of precious observing time. In- 
deed, we may have viewed such a scenario in one 
of our stars GJ1253. 

Fig. [7] shows two epochs of spectroscopic mea- 
surements of the region around the Ha line for the 
star GJ1253. The upper panel is the first obser- 
vation, made on 08/30/2008, whereas the lower 
panel shows the same region only one month later 
(09/26/2008). In the first observation there is no 
indication of any Ha emission present, however 
only 27 days later when the second measurement 
was obtained, there is clear evidence for the Ha 
line in emission. Such a short period evolution of 
the line, coupled with the relatively low v sin i of 
the star (< 4.5 km s _1 ) may indicate that a small 
magnetic region is rotating in and out of view. If 
confirmed continuous monitoring of such Ha line 
changes could lead to the true rotation period of 
the star, removing the sin i degeneracy and giv- 
ing both the true velocity and, for any discovered 
companions, their true masses without the need 
for continuous pho tometric monitoring. Note that 
Byrne et al.l (|l996l ) found variable Ha emission in 



the active M star HK Aquarii which they attribute 
to prominence like clouds above the stellar sur- 
face that have characteristic timescales less than 
half the photometrically determined rotation pe- 
riod. It may be the case that such variability of Ha 
flux is mainly generated through flares or a chang- 
ing magnetic field structure. Indeed, Ha vari- 
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ability for M stars is significant across the spec- 
tral domain, ev en on rela t ively short time scales 
(<60 minutes). iLee et al.l (|2009l ) find that ~80% 
of their mid-to-late M star sample exhibited sta- 
tistically significant Ha variability. Therefore, if 
any correlation does exhist between Ha variabil- 
ity and rotation velocity the characteristics must 
be sufficiently large such that typical short period 
small scale variations do not mask it out. 

Finally, we also searched for significant Ha ab- 
sorption and found eight stars (14%) that show 
significant absorption. All of these have v sin is 
that would meet our criteria for planet search se- 
lection, as would be expected from a rotation- 
activity connection. Only, 2 out of the 8 objects 
that exhibit significant Ha absorption are below 
our resolution limit, with a further 3 in agreement 
with the adopted resolution limit to within the 
estimated total uncertainty of the analysis tech- 
nique. All stars with no Ha emission would meet 
our criteria for planet search selection, with the 
largest rotation rate of these coming from the M5 
dwarf GJII82 (v sin i=6.8±0.5 km s" 1 ). Col- 
umn f 3 of Table [3] shows the Ha emission flag, 
where either a detection was made or not, and 
stars with no significant emission which have as- 
terisks show evidence for Ha absorption. 

3.3. Spectroscopic Binaries 

There is the possibility that stars with miss 
shapen profiles are spectroscopic binaries of simi- 
lar spectral type or stars with powerful magnetic 
fields. Fig. [8] shows blended profiles found for both 
GJf 080 (upper panel) and LHS3080 (lower panel), 
ft is clear that for both of these systems there is 
a primary deconvolved profile and a weaker sec- 
ondary. To test if these are real rather than de- 
convolution artifacts, we split the spectra into two 
wavelength ranges, one spectrum covering the blue 
orders and one covering the red orders. By decon- 
volving both the red and blue spectra indepen- 
dently, we were able to confirm that these stars 
still exhibit double line profiles, lending weight to 
fact that the double profiles are due to real phe- 
nomena. The spectra for GJ1080 shows double- 
lined profiles for atomic and molecular lines. The 
double-lined profiles for LHS3080 were difficult to 
confirm by eye in the spectrum. Along with both 
of these, there are another two binary systems, 
those of GJ3f29 and G1802. Both of these stars 



have double profiles, but with much wider sepa- 
rations, even though they exhibit noisier spectra. 
These are more clear cases than the previous two 
since they also have double Ha emission profiles. 
G1802 was known to have a low-mass com- 



(Pravdo et al. 


20051). 


by Llovd et al. 


(20061) 



which was di rectly imaged 



In addition, llreland et al 



(|2008l ) has recently found this star to be part of 
a triple system, with a short period (P~19 hrs) 
spectroscopic companion. Our final deconvolu- 
tion for this star is very noisy, in part due to 
the blended light from three different components 
with different spectral profiles. We do note how- 
ever that there are two strong profiles in the fi- 
nal deconvolution (properties listed in Table [T|) 
and potentially a further two weaker profiles. The 
weak profiles are only borderline significant due 
to the associated noise, but one might explain the 
short period binary found by Ireland et al. and the 
other, which is widely separated from the profile 
of G1802A, could be an additional longer period 
companion. 

McCarthy fc Zuckermanl (|2004l ) surveyed both 



GJ3I29 and LHS3080 as part of their AO cam- 
paign to detect substellar companions to a host 
of nearby young M dwarfs. They found no vi- 
able candidate to these stars between around 5" 
to 15" arcseconds. The lack of any viable detec- 
tion around these, apparently young, stars is prob- 
ably explained by the close separation of the ob- 
jects, since the detectability of spectroscopic bi- 
naries are skewed towards short period compan- 
ions. Indeed, the LHS3080 profiles are blended 
with each other and since the secondary profile is 
significantly weaker than the primary, it is prob- 
ably a much later and fainter companion object. 
The properties of these spectroscopic binaries are 
listed in Table [f] which include the Hipparcos V, 
2MASS J, H and K s apparent magnitudes, their 
spectral types, the velocity separation of the dou- 
ble peaks along with their combined uncertain- 
ties, the v sin i estimates for both the primary 
star (v sin i pr ;) and the secondary (v sin i sec ) and 
the Ha flag. Note the small uncertainties are an 
artifact of the blended fitting procedure and do 
not necessarily mean these values are more pre- 
cise than the single profiles. 
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Table 1 : M dwarf binary candidates 



Star 


V 


J 


H 


K a 


Spec Type 


Av (km s~ x ) 


v sin ip r i (km s 1 ) 


v sin i BCC (km s L ) 


Ha Emission 


GJ1080 


12.81 


8.98 


8.50 


8.22 


M3.0 


28.44±0.28 


<3.0 


<3.0 


No 


GJ3129 


14.27 


9.65 


9.06 


8.80 


M4.5 


89.36±0.78 


5.2±0.1 


6.6±0.3 


Yes - Double 


G1802 


14.67 


9.56 


9.06 


8.75 


M5.0 


150.03±1.19 


6.4±0.4 


<4.5 


Yes - Double 


LHS3080 


14.28 


9.67 


9.11 


8.82 


M4.5 


18.35±0.32 


<3.0 


3.8±0.1 


Yes - Single 
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Fig. 7. — Region around the Ha line profile for the 
star GJ1253 at two separate epochs. The top plot 
shows the first measurement made on 08/30/2008 
and no Ha emission was found. However, the 
lower plot shows the second measurement made 
for this star one month later (09/26/2008) and 
clear emission is present. 
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Fig. 8. — Example blended deconvolved profiles 
for two of the probable binary systems in this sam- 
ple. The upper plot represents the star GJ1080 
and has two clearly defined deep profiles, whereas 
the lower plot, which represents the LHS3080 sys- 
tem, has a secondary profile also, possibly from a 
weak signal contribution from the secondary. The 
dashed curves represent the best fit double profiles 
to the data. The two other potential binaries have 
profiles with much larger separations. 
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4. Temperature, Mass and Metallicity 

When performing radial-velocity searches for 
exoplanets other characteristics of the star deter- 
mine radial-velocity limits e.g. the stellar mass. 
The mass of the primary determines the lower 
limit to the planetary mass for a given radial- 
velocity amplitude and therefore knowledge of 
the mass of any planet search target star is es- 
sential. We have determined the stellar mass 
for the bulk of our sample by utilising relations 
between M star mass and their absolute mag- 
nitudes. To generate accurate absolute magni- 
tudes we se arched the Yale Trigon ometric Paral- 
lax project (Ivan Altena et ai1 ll995l). the RECONS 
list ([Henry et al.ll2006l ) and any parallaxes in the 
Gliese Catalogue of Nearby Stars to obtain ac- 
curate parallax measurements for all our candi- 
dates. The parallaxes, along with their associ- 
ated uncertainties, are shown in column 8 of Ta- 
ble [3] and these were used to determine the abso- 
lute magnitudes by converting them to distance 
and measuring the distance modulus. The pho- 
tometry was acquired using both the SimbaoH 
and Viziejf) astronomical databases. The near in- 
frared p hotometry was taken f rom the 2MASS cat- 
alogue (jSkrutskie et al. 120061) and the K s mag- 
nitudes were converted t o K using the calibra- 
tions in ICarpenteri (|200lh . The absolute V and 
K magnitudes were t hen input into the em piri- 
cal mass relations from Dclfossc ct al. (2000), and 
two measurements of the mass for each star were 
determined. Both these measurements were av- 
eraged to get the final combined mass estimates 
for all objects with known parallax and their un- 
certainties were taken as the standard deviation 
of the two measurements. The uncertainties are 
typically ±10-15%, with a few of the closest stars 
having uncertainties down to the ±1% level. Both 
the masses and their uncertainties are shown in 
column 9 of Table [D 

Along with the mass, both the effective temper- 
atures and metallicities ([Fe/H]) were estimated 
for the sample. The effective temperatures were 
determined using t he V — K s relation taken from 
Casagrande et al.l (|2008l ) and has a typical inter- 
nal uncertainty of ±17K, however the uncertainty 
on the overall accuracy of the technique will prob- 



ably be substantially larger than this. The val- 
ues with question marks next to them flag highly 
suspect effective temperatures. The [Fe/H] abun- 
dances were firstly determined photo metrically us 



i ng th e colour magnitude relation in lBonfils et al 
(|2005l ). They quote the typical uncertainty for 



this method as ±0.2 dex which we assign to all 
our metallicities shown in column 10 of Table [3] 
Casagrande et al. also provide metallicities using 
their method and these agree with the Bonfils et 
al. values wi thin the quoted uncerta inties. Note 
that recently IJohnson fc Appd (|2009h claim that 
the Bonfils et al. relation may underestimate the 
metallicities of M dwarfs by as much as 0.3 dex. 
Finally, Table [3] also shows each star's V, J, H and 
K s photometry. 

5. Results: v sin i Distributions 

Fig. [5] shows the distribution of early-to- 
late M star rotation rates against their spec- 
tral type. The filled circles represent the stars 
in this work and the o pen circles represent 
litera t ure v sin is (from IStauffer fe Hartmann 
19861: iMarcv fe Chenl fl992t iDelfosse et alJ Il99; 



Gizis et all200a iMohantv fe Basrl2 003: Ba iler-Joned 
2004 iFuhrmeister fc Schmittl 12004 Ijones et al. 



http://simbad.u-strasbg.fr/simbad/ 



2 Simbad website: 
3 Vizier website: http://webviz.u-strasbg.fr/viz-bin/VizieR 



2003: lReinersll2007t IWest fc Basrill2009h . All lit 
erature v sin is are shown in Table 0] along with 
any measured V magnitudes, 2MASS J, H and 
K s photometry, spectral types, effective tempera- 
tures, parallaxes, masses and [Fe/H] abundances, 
all determined using the same methods outlined 
above. 

There is a large spread in v sin i across the mid 
M star regime (M4.0-M6.5) compared with early 
M stars (<M4). The mid M objects range from al- 
most as high as over 50 km s _1 and down as low as 
essentially zero km s _1 . The median dispersion of 
the v sin i distributions in the spectral type bins 
between M0-M3 is 3.70±0.79 km s _1 , compared 
to the median dispersion of 9.00±8.31 km s _1 for 
bins at M4-M9.5. However, if we allow for a spec- 
tral typing uncertainty of ±0.5 sub-types and ex- 
clude both the M0 and M9.5 spectral bins, since 
they could be contaminated by K and L dwarfs, 
we find the dispersions are 3.70±0.87 km s _1 and 
9.00±4.37 km s _1 respectively. A K-S test is run 
to compare these distributions and this returns a 
D-statistic of 0.539, giving an extremely low prob- 
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Spectral Type (M) 

Fig. 9. — The distribution of v sin is against mea- 
sured spectral type plotted in log step. The val- 
ues from this data are shown by filled circles, with 
the literature values shown by open circles. The 
trend of increasing rotation velocity with decreas- 
ing temperature is seen here and is highlighted by 
the solid line linking the filled stars. These stars 
mark the median values for all data in each spec- 
tral bin. The uncertainties plotted represent Pois- 
son statistics and the downward pointing arrows 
mark upper limits. Represented by the dashed 
lines are the best straight line fits to the medians 
for the early M dwarfs and mid-to-late M dwarfs, 
along with their standard ±lcr uncertainties plot- 
ted by the dotted lines. The early Ms are found to 
have a much flatter distribution than the mid-to- 
late Ms, which are tending towards a more rising 
trend. 



ability of only 1.396xl0 _8 % that the rotation rates 
for stars between M0.5-M3 are drawn from the 
same parent population as those between M4-M9. 
Also, a large fraction of stars with spectral types 
below M6.5, and particularly below the convec- 
tive boundary region, only have measured upper 
limits. This will bias this result towards a non- 
correlation and so we can expect that with fur- 
ther detections at low velocities for these stars, 
this correlation will become more pronounced. 

This step from low to high dispersions be- 
tween M3 to M4 spectral types is thought to 
be due to the increased spin-down timescale 



towards decreasing m ass ( Delfosse et al. 1998; 



Mohantv fc Basril 120031 ) and that the spin-down 
timescale is a significan t fraction of the a ge of 
the young d i sk. In deed, iDonati et al (2008) and 
Morin et al. ( 2008h have shown that the magnetic 



field topologies of early M's (M0-M3) and mid- 
M's (M4) are significantly different, with the early 
M's having mainly toroidal and non-axisymmetric 
poloidal fields, whereas the M4s mainly exhibit 
axisymmetric poloidal fields. This result indicates 
a change in the magnetic field properties at the 
classical boundary between partially radiative and 
fully convective envelopes. If the spin down times 
are governed by the magnetic fields in this regime 
then a differing v sin i distribution might be ex- 
pected. It might be the case that axisymmetric 
poloidal fields interact with the stellar wind more 
weakly than the toroidal, non-axisymmetric fields, 
driving a less efficient braking mechanism. Inter- 
estingly, a probe of the spindown timescale and 
the mechanism driving it can be made if we take 
the Delfosse et al. and Mohanty & Basri v sin i 
and log(LH cl /LBoi) values (Mohanty & Basri show 
the activities in these two works to be in excellent 
agreement). We can use these values to trace the 
rotational history for the fully convective, mid- 
M star regime (M4-M7). When the stars are 
binned into their respective spectral types, the 
trend found between these two quantaties reveals 
a sharp saturation boundary (e.g. Fi g. 9 in Mo- 
hanty & Basri). By using Eq 11 . 3.1 in IWest et al.1 



(2009), one can use the activity values to gain a 
statistical insight into the age distribution of these 
samples. In particular, the saturation boundary 
changes as a function of spectral type, increasing 
in velocity with increasing spectral type. Table [2] 
shows the spectral type bins, along with the sat- 
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uration boundary in velocity and activity lifetime 
(I) given in West et al. Due to the spindown 
of stars by the wind braking mechanism we can 
use this table to say that from a statistical foot- 
ing, the average age of M4 stars with v sin is 

< 8 km s _1 is older than 5 Gyrs, M5s with ve- 
locities < 10 km s -1 , and M6s with velocities 

< 12 km s -1 , are older than 7 Gyrs and M7s with 
velocities < 13 km s _1 are older than 8 Gyrs. 

Fig. [TU] shows the v sin i values against spec- 
tral type for all stars in this study with a split 
made by their Ha status. The filled circles rep- 
resent all stars with no significant Ha emission or 
absorption, the open circles represent stars with 
Ha absorption and the open stars represent ob- 
jects with Ha emission. The rotation-activity con- 
nection suggests a correlation between Ha emis- 
sion and the rotation velocity of cool M stars and 
this is clearly seen in this plot. All stars with 
v sin i > 7 km s~ x are found to exhibit signif- 
icant Ha emission. Indeed, we have just shown 
that, on average, mid-to-late M stars with veloc- 
ities > 8 km s _1 are still in their young, active 
phase of evolution. No stars later than M5 were 
found to exhibit Ha absorption, with five of these 
exhibiting no significant Ha emission (55%). Such 
numbers correlate with previous results which sug- 
gest a high frequency of act ive M stars towards 
the latest spectral ty pes (e.g. Fleming et al. 200(A 
iMohantv et al.ll2002h . even though we find a simi- 
lar fraction of active M stars in all of our spectral 
bins. However, in our final M6.5 bin both stars are 
found to exhibit Ha emission which agrees with 
the increasing trend of activity towards values ap- 
proaching 100% at a spectral type of around M7. 
This increase in stars with Ha in emission towards 
later spectral types appears to be a product of 
the fraction of young disk stars in this regime and 
also the increase in activity lifetime with increas- 
ing spectral type. From Table [2] we see that both 
the activity lifetime and the time spent in a state 
of high rotation increases with increasing spectral 



Table 2: M dwarf activity saturation limits 



Spectral Type 


v sin i (km s 1 ) 


I (Gyr) 


M4 


8 


5 


M5 


10 


7 


M6 


12 


7 


M7 


13 


8 



type. Such trends would naturally give rise to a 
larger fraction of stars with Ha profiles in emission 
towards the latest M stars. 

The filled stars in Fig. [9] represent the median 
values in each spectral bin along with their asso- 
ciated Poisson errors. The solid line connects the 
points and visually highlights t he increasin g trend 

(2002) 



Reid et al 



towards later spectral types 
suggest a flat distribution of rotation rates be- 
tween stars of M6-M9, whereas larger samples in- 
dicate this may not b e the case (Mohanty et al.; 
iReiners fc Basril 2008). The medians indicate a 
rising trend through the mid-to-late M dwarfs, 
flattening off towards the end. Due to the change 
in distributions at the fully convective boundary 
we perform two straight line fits to the medians 
in the spectral ranges M0-M3.5 and M4-M9. The 
best fit to the early M stars clearly shows a flat 
trend across the whole regime, whereas the fit 
to the late Ms show a rising trend, appearing to 
flatten towards the later M stars (note the cur- 
vature due to logarithmic plotting). The fits to 
the early and late M star samples are described 
by v sin i=0.09(±0.30)xSpT±3.41(±0.67) km s" 1 
and v sin i=2.10(±0.53)xSpT-4.54(±3.56) km s"\ 
with standard uncertainties, represented by the 
dotted lines in the plot, of ±0.89 km s -1 and 
±2.79 km s _1 respectively. The total sample of 
less than 300 is still rather small, especially when 
binned. Indeed there appears a dearth of ob- 
jects between M6.5-M8.5 with rotation rates above 
~15 km s , which may indicate another popula- 
tion change above M6.5. The evidence for this 
gap is weak at present due to the low number of 
stars in these spectral bins, therefore further ob- 
servations and are needed and any biases studied 
to fully validate the existence of this feature. 

At temperatures below around 2800K (ap- 
proximately M6-M7 type objects) dust formation 
and opacity are important in st ellar/substellar 
atmo s pheres (Tsuii et al.l Il996t I Jones fc Tsuii 



1997 



199a IChabrier et aiT[ 2000: 



Tinney et al 

Baraffe et al.ll2002 and references therein) . iBerger et al 



(2008) have shown that late-M stars mark a transi- 
tion in the properties of the magnetic field and its 
dissipation, along with high temperature plasma 
being generated in the outer atmosphere. They 
go on to hypothesize that the stellar rotation may 
play a part in this process and indeed the differ- 
ence shown here between the mid and late type 
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M stars seems to add to this conclusion. A K-S 
test reveals a D-statistic of 0.639, or 5.413xl0~ 6 %, 
that stars in the range M0.0-M6.5 and those in the 
range M7.0-M9.5 are drawn from the same parent 
distribution. However, given that we have already 
shown there to be a large difference between early 
Ms (M0-M3.5), this will bias this probability test. 
When we remove all stars earlier than M4 we find 
a D statistic of 0.544, which relates to a probabil- 
ity of only 1.437xl0 _3 % that these are drawn from 
the same parent population. This >5<t result may 
indicate that at temperatures when dust opacity 
becomes important there is a change in the ro- 
tational braking mechanisms and hence the mag- 
netic properties of ultracool dwarfs. This might 
give rise to the flattening trend indicated between 
M6.5-M9 stars, however a more comprehensive 
study is needed, particularly to decouple the age 
of these stars by studying the space motion to de- 
termine if they are young or old disk stars. Also 
the biases of the literature surveys are important. 
For in stance, studies like those of IWest fc Basri 
(2009) 1'ocus only on selecting inactive, and hence 
slowly rotating, late type M stars. In addition, 
current models show that the late M star regime 
can also be populated by young brown dwarfs. 
Finally, this relation also suffers from the lack of 
low v sin i detections already mentioned above, 
even more so given the M6.5 detection bound- 
ary. Therefore we expect this result might also 
become more pronounced with further low v sin i 
detections at spectral types below M6.5. 

The normalized distribution of v sin is are rep- 
resented by the histograms in Fig. QTJ where the 
solid histogram is for all stars in the spectral range 
between M0-M3.5 and the dashed histogram is 
for all stars in the range M4-M9.5. These in- 
clude all values determined in this work com- 
bined with those in the literature. It is appar- 
ent that both distributions peak at low rotation 
rates (~3 km s _1 ), with peak values of 55 and 42 
stars respectively. We find that the total num- 
ber of v sin is < 10 km s _1 is 198 and these 
should represent useful stars for future near in- 
frared rad ial-velocity pl a net se arch projects such 
as PRVS. iBouchv et all (|200lh show that the in- 
formation content drops by a factor of ^3.5 be- 
tween rotation velocities of ^2-10 km s -1 , mak- 
ing > 10 km s _1 a reasonable M star radial- 
velocity selection cut. This sample is still large 



(124) when we include all mid-to-late M stars in 
the range M3-M9.5 (stars where obtaining optical 
precision radial-velocities becomes extremely diffi- 
cult). Note the binary systems have been left out 
of Figs. [TT] and O since the combined luminosities 
will generate inaccurate photometry and therefore 
inaccurate spectral types. Also, binary systems 
like these make radial-velocity exoplanet searches 
much harder since any small planetary signature is 
masked by the large short period binary velocity, 
meaning these are not ideal planet search targets 
for precision radial- velocity programs moving into 
an unexplored parameter space. 

Comparing the distributions of both histograms 
helps to probe the possible changing rotation prop- 
erties of M stars at the fully convective boundary. 
We have employed two power law fits to each dis- 
tribution separately in order to test the changing 
velocity distribution between these two regimes. 
The red (dark grey)0 curve is fit to the sample 
of early M dwarfs between M0-M3.5, whereas the 
green (light grey) curve is fit to the mid-to-late M 
dwarf sample. The fits are made to the bins by in- 
cluding Poisson uncertainties which are not shown 
in the plot for clarity. It can be seen that the fit to 
the early M stars drops much more rapidly than 
the fit to the later Ms. The power laws are de- 
scribed by SN/Sv sin i oc x~ 313 for the early Ms, 
whereas for the later Ms it is only found to be 
oc x -112 , highlighting the differing steepness of 
each slope. The early Ms have a much longer tail 
than the later Ms due to this faster decay of the 
distribution and this is also probably an underes- 
timate since we included all upper limits in the 
data to increase the sampl e size, which included 
a num ber of stars from the IStauffer fc Hartmannl 
( 19861 ) sample with detection limits of 10 km s _1 . 
We note that the curvature of these slopes are also 
affected by the activity lifetimes shown in Table [2] 
since they change with spectral type. 

The difference of 2.01 in the exponent between 
the two spectral samples allows us an insight into 
the efficiency of the braking mechanism between 
partially and fully convective stars, assuming the 
larger measured rotation in the later M star sam- 
ple is not due to increased line blending from 



4 The colours in brackets relate to the printed document, 
whereas the colours in the text are for the online edition of 
the article. 
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increased molecular bands and instrument al res- 
olution arguments. iDelfosse et al have 
shown that the majority of rapidly rotating mid-M 
stars are members of the young disk population, 
whereas the older population tend to rotate more 
slowly. If this is indeed the case, then the brak- 
ing mechanism is at play in later M stars but the 
efficiency has dropped across the fully convective 
boundary. As mentioned earlier, the change in 
the field topology between the partially and fully 
convective boundary is probably the driving fac- 
tor which governs th e efficiency of the wind brak- 
ing mechanism since iReiners fc Basri ( 2007 ) have 
shown that fully convective stars produce field 
strengths as strong as partially convective stars. 
To better probe the braking mechanism in this 
fashion, in addition to gaining more data, it is nec- 
essary to also fold in an age proxy for the sample, 
and decouple both the young and old disk popula- 
tions to compare these objects. Along with this a 
better understanding of the activity lifetimes and 
how these change with spectral type should be 
considered. More magnetic field topology stud- 
ies are required for later M stars which can help 
to confirm if the later Ms also have axisymetric 
poloidal fields and add weight to the topology ar- 
gument. Finally, more detailed testing of changes 
around the dusty regime (~M6.5) where stars can 
also be young brown dwarfs might usefully be in- 
vestigated. 

6. Conclusions 

We present the initial results from our study of 
rotation rates for a range of M stars as part of 
our target selection for a PRVS-like planet search 
project. We observed over 50 M stars with HRS 
on the HET with the aim of selecting the slow- 
est rotators in order that a near infrared planet 
search survey, such as PRVS, shall have a statisti- 
cally large sample of M stars where highly precise 
radial-velocity measurements can be accrued. Of 
our sample of 49 suspected single M stars between 
M3-M6.5, we find 36 have v sin is < 10 km s _1 
which will represent good radial-velocity targets. 
When we include all literature M stars in the opti- 
cally difficult radial- velocity regime (M3-M9.5) we 
find this value increases to 124. 

We also confirm the increase of rotational veloc- 
ities between early to mid M stars. This change at 



the fully convective boundary seems to be linked to 
a change in the topology of the magnetic fields be- 
tween such stars, indicating axisymetric poloidal 
fields drive a less efficient wind braking mecha- 
nism. Also there is an initial indication that stars 
with spectral types in the range M6.5-M8.5 have a 
different rotational velocity distribution compared 
with those of below M6.5, as the distribution ap- 
pears to flatten off beyond this regime. Since this 
is around the temperature where dust opacity be- 
comes important (^2800K), there may be another 
change in the efficiency of the braking mechanism 
in such stars which could indicate another mag- 
netic field topology change. In addition, we also 
show how knowledge of the v sin i can be used 
to put a lower, or upper, limit on the age of mid- 
M dwarf stars, since the rotation-activity relation 
has a temperature dependent saturation level. 

We also highlight the rotation-activity relation 
through emission, or lack thereof, of the Ha line. 
There appears a boundary of around 7 km s _1 be- 
tween stars with and without Ha emission, with 
the fast rotators almost always exhibiting such 
emission, however since the sin i degeneracy is 
present it is difficult to account for any firm empir- 
ical boundaries with small numbers. We observed 
the star GJ1253 over two epochs and found the 
Ha emission had switched on over a period of less 
than one month. This may be due to an active 
region rotating in and out of our field of view with 
a period equal to the rotation period of the star, 
or due to a flaring event. 

We have also discovered three spectroscopic bi- 
nary systems and confirmed another. Both com- 
ponents in the GJ1080 system produce fairly sim- 
ilar profiles, indicating they are of similar spectral 
type, the secondary likely a little cooler and less 
luminous. The profiles for GJ3129 and G1802 are 
widely separated and both of these exhibit dou- 
ble Ha emission features. G1802 is also a known 
triple system, but we find evidence for another 
companion in the system, however given the noise 
due to blended light from at least three separate 
sources the evidence is weak. In comparison to 
these fairly strong profiles the secondary profile in 
the LHS3080 system is significantly weaker, indi- 
cating the companion is significantly cooler than 
its host star. Finally, we have flagged other M 
dwarfs where there is some evidence for binary 
companions in these systems. 
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Table 3: Characteristics of M stars in this study. 



Cn 



Star 


V 


.7 


H 


K s 


Spec Type 


t E ff ( k ) 


tt (ma,) 




M/M ( 




[Fe/H] phot 


V sin « (km s- 1 ) 


FWHM M (km s 1 ) 


Ha Emissioi 


G121-028 


14 




9.94 


9.31 


9.02 






2948 


54.00±8.00 





200±0 


005 


- 


3.8±0.7 


7.2 


No* 


G180-011f 


13 
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The f symbol relates to deconvolved profiles that show cither weak evidence for a possible binary component, but not significant enough to be included in the binary table, however it may affect the 
deconvolved profile enough to give rise to a larger rotation velocity than the star has. or stars where the dcconvoluted profile has low S/N and the construction is of lower quality, particularly in the 
wings, which again can give rise to an inaccurate value for the v sin i. The question marks after the effective temperature arc used to flag highly suspect values. The asterisks in the Ha column represent 
stars with significant absorption. The uncertainties shown for the v sin is represent the formal uncertainties on each parameter fit. However, as explained in the text, the actual uncertainties for high 
S/N spectra of objects rotating <20 km s~~ ^ is zb 1 km s ^ , increasing to around il0-20% for the fastest rotating objects. 



Table 4: Tabic of v sin is taken from the literature and used in this work. Columns are the same as Table [3] 



except without the telluric widths in column 15. 
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Fig. 10. — The distribution of rotational velocities 
against spectral type for stars in this study with 
a split based on each star's Ha profile. The filled 
circles represent stars with no significant Ha emis- 
sion or absorption. The open circles are stars with 
significant Ha absorption and the stars represent 
objects with detected Ha emission. All stars with 
velocities greater than 7 km s^ 1 have detected Ha 
emission. Also no significant Ha absorption was 
detected in stars later than M5. 




O 10 20 30 40 50 60 

v sin i (km/s) 

Fig. 11. — Histograms of rotation velocities in 
this sample and in the literature split by spectral 
type. The solid curve represents the stars in the 
spectral range from M0-M3.5, whereas the dashed 
curve represents the stars between M4-M9.5. Both 
samples peak at low rotation velocities of around 
^3 km s _1 , however the bins that contain the later 
type objects have many more stars with measure- 
able v sin i. The solid curves are the best fit 
power laws to the data, with the red (dark grey) 
curve representing the M0-M3.5 bins and the green 
(light grey) curve representing the M4-M9.5 data. 
The changing power law between the two spec- 
tral regions highlight a possible change in the ro- 
tational distribution for fully convective stars. 
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